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Detection of CO in the inner part of M31's bulge 
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ABSTRACT 

We report the first detection of CO in M31's bulge. The 12 C0 (1-0) and (2-1) lines 
are both detected in the dust complex D395A/393/384, at 1.'3 (~ 0.35 kpc) from the 
centre. From these data and from visual extinction data, we derive a CO-luminosity 
to reddening ratio (and a CO-luminosity to H2 column density ratio) quite similar to 
that observed in the local Galactic clouds. The (2-1) to (1-0) line intensity ratio points 
to a CO rotational temperature and a gas kinetic temperature > 10 K. The molecular 
mass of the complex, inside a 25" (100 pc) region, is 1.5 10 4 Mq. 
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1 INTRODUCTION 

The bulk of gas of M31 lies between 6 and 18 kpc and fol- 
lows a pattern of thin spiral arms (e.g. Dame et al. 1993, 
Neininger et al. 1998) . Despite the presence of dust lanes, 
massive stars, and evolved stars (e.g. O'Connell et al. 1992; 
Davidge et al. 1997), there is little HI atomic gas in the 
inner bulge and, up to now, molecular gas had escaped de- 
tection: Sofue & Yoshida (1993) reported the detection the 
12 CO (1-0) line in the dust complex D395/393 at less than 
1 kpc from the centre, but this detection was invalidated by 
Loinard, Allen & Lequeux (1996) who re-observed the same 
complex with a much better sensitivity. So, the detection of 
CO in the dust complexes D 478 and D 268, at 2-4 kpc from 
the central by Allen & Lequeux (1993) and Loinard & Allen 
(1998), was the closest to the centre ever reported. Noting 
that the CO (1-0) to (2-1) intensity ratio was lower in these 
clouds than in the Galactic Giant Molecular Clouds, these 
authors concluded that CO may be so cold in the inner disc 
of M31, that its excitation temperature barely exceeds the 
cosmic background temperature. They thus suggested that 
although molecular gas and CO could be abundant, the CO 
mm lines are so weak that they escape detection. 

In order to better understand the gas budget at the 
centre of M31, we have embarked on a reanalysis of the op- 
tical data. We found that the positions observed by Loinard, 
Allen & Lequeux (1996) were not centred on the darkest dust 
patches, which prompted us, taking advantage of the instal- 



* We assume a distance to M31 of 780 kpc (e.g. Holland 1998; 
Stanek & Garnavich 1998), i.e. l'.'=3.8pc. 



lation of new generation mm-wave receivers on the IRAM 
30-m telescope, to make new CO observations. 



2 OBSERVATIONS 

2.1 Distribution of the extinction in the bulge 

We used a compilation of optical images of M 31 with res- 
olutions (FWHM) close to 1 arc-sec to locate regions with 
significant absorption in the bulge of M31: (1) a B frame 
(O'.'9/pixel) covering a field of 22f9 x 17.'9 obtained during 
the course of a 5 yr nova search by Ciardullo et al. (1988); (2) 
B and R images (O'.'3/pixel), covering a field of 12' x 9', col- 
lected during an intensive monitoring of M 31's bulge aimed 
at micro- lensing detections (Ansari et al. 1997, 1999), here- 
after the AGAPE data; (3) BVRI images from J. Tonry 
covering a 4' x 4' field (Metzger, Tonry & Luppino 1993). 
The astrometry (M. Auriere, private communication) of the 
optical data presents small (~ 1 arcsec) distortions on the 
degree scale; this is insignificant given the size of the beam of 
the radio telescope (21" at 115 GHz and 11" at 230 GHz). 
We model the photometry of M31's bulge with elliptical 
annuli using the standard surface photometry algorithm de- 
veloped for IRAF (Jedrzejewski 1987). This model intends 
to reproduce the light profile along M31's bulge without 
extinction. Hence, the ellipse geometry parameters (centre, 
ellipticity, position angle) are fitted in this procedure. In 
addition, to avoid possible perturbations of the fit due to 
extincted areas, only the median intensity over the ellipti- 
cal annulus sectors is used. The observed extinction is then 
defined as Aa = — 2.5 log lo (0obs/0modei) where </> bs is the 
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Figure 1. The Ab extinction obtained from the data of Ciardullo et al. (1988) is shown in grey-scale. The black contour levels correspond 
to 0.1, 0.2 and 0.3 mag in Ag scale. The two straight segments intersect at the centre of M31. Their lengths correspond to 1 kpc along 
the major and minor axis. The circle shows the ON position toward the dust complex D395A(/393/384). Its diameter is the telescope 
HPBW at 115 GHz. The small crosses display the OFF positions used during the observations. The red contour levels give the possible 
contamination level in the OFF position at 115 GHz, assuming CO emission is proportional to visual extinction (first 0.05 per cent, 
second 0.10 per cent, third 0.15 per cent). 



observed brightness and <^ mo d e i the brightness in the model. 
This derivation of Aa, which assumes that the dusty cloud 
is in front of the bulge, provides a lower bound of the true 
internal extinction. 

The comparison of 3 completely independent B data 
sets shows that no significant observational artifacts are in- 
troduced in the extinction maps. The nominal centre of M 31 
is taken as the optical nucleus at aj2ooo.o = 00 h 42 m 44 s . 371 



and 5j2ooo.o = 41° 16' 08734 following Crane, Dickel & 
Cowan (1992). Figure shows the extinction obtained in B, 
centred on the dust patch D395 A/393. All the structures 
seen here are well correlated with those seen in B — R. We 
note that the procedure defining the ellipses tends to remove 
any smooth structures within 10 arcsec from the centre. The 
structure in the north-western part is the edge of a much 
larger extinction pattern; the most prominent structure de- 
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Figure 2. The telescope HPBW at 115 GHz and 230 GHz during 
ON scans (a J2 000.0 = 00 h 42 m 42M0 and <5j 20 00.0 = 41°17'24'.'l), 
superimposed on the Ag extinction map derived from the 
AGAPE data. The small circle on the left marks the position 
observed by Loinard et al. (1996). 

tected near the centre is D395 A/393/384 with a typical size 
of 95 pc (25 arcsec). It is centred at —26", +76" from the 
nucleus, and consists of several dark patches or clumps. 

2.2 Radio observations 

The observations were carried out on 1999 June 13-14 with 
the IRAM 30-m telescope. Most of the observing was made 
in the symmetrical wobbler switching mode where the sec- 
ondary mirror nutates up to a maximum limit of ±240" 
in azimuth. The beam throw was determined as a function 
of the hour angle in such a way that OFF positions lie in 
extinction-free regions (see Fig. 111). Near transit, we had to 
use position switching mode, taking an extinction-free OFF 
position located at 122", 321" from the nucleus. Figure^ 
show the areas covered by the telescope beam (HPBW) in 
the ON position at 115 and 230 GHz. The beam is centred 
as best as possible on the dust complex, and encompasses 
most of D395A ( and also, at 115 GHz, parts of D393 and 
D384). The telescope pointing was checked every hour on 
nearby quasars and found to be accurate within 3 arcsec. 

We used four receivers simultaneously, two for 12 CO(l- 
0) at 115 GHz and two for 12 CO(2-l) at 230 GHz. Each 
115 GHz receiver was connected to two autocorrelator sub- 
bands, shifted by 40 MHz from each other. Each sub-band 
consisted of 225 channels separated by 1.25 MHz. A very 
simple algorithm, combining such a pair of settings, was 
used to remove occasional "platforming" between the 80 
MHz wide units which compose the sub-bands. Each 230 
GHz receiver was connected to a filter-bank consisting of 
512 1-MHz-wide channels. Linear baselines were fitted and 
subtracted from the spectra registered every 60 s, after 
the platforming correction for the data at 115 GHz. Some 
spectra were flagged, especially for one of the receivers at 
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Figure 3. 12 CO(1-0) and 12 CO(2-l) spectra at a resolution of 3.0 
km s _1 , observed towards the dust complex D395A(/393/384) 
(offsets from M31's nucleus: —26", 76"). The r.m.s. noise levels 
in individual channels are illustrated by the dashed lines; they 
correspond respectively to 1.8 and 2.9 mK. The integration times 
(IT) quoted correspond to the sum of the ON and OFF times 
for two receivers. IT is shorter for the (1-0) spectrum, because 
we excluded the data with OFF positions closer than 2 HPBW 
(50") from the ON position. The equivalent system temperature 
is calculated from the r.m.s. noise and IT, using the standard 
radiometer formula. 

230 GHz, due to temporary instabilities. The spectra were 
then combined, using inverse variance weights, which corre- 
spond to the residuals about the baselines. 

The resulting spectra are displayed in Fig.^j and the 
lines parameters are summarised in Table [j]. We estimate 
that even if the gas in the OFF positions had the same ve- 
locity as in the ON position and if the CO luminosity was 
proportional to A_b, we would have decreased at most by 6 
per cent the 115 GHz line intensity (see the red contours in 
Fig.hL The uncertainties on the line intensities quoted in 
Tablehl account for the r.m.s. noise and for baseline uncer- 
tainties; they do not include calibration uncertainties which, 
according to a check on the standard source DR21, should 
be < 10 per cent. (We adopted a T^ flux scale corresponding 
to integrated 12 CO(1-0) and 12 CO(2-l) emissions of respec- 
tively 410 K km s" 1 and 360 K km s" 1 at the centre of DR21 
- see Mauersberger et al. 1989). We note that, although 
the spectra of Fig. H were observed only 15" away from those 
of Loinard, Allen & Lequeux (1996), the line intensities are 
much weaker and the velocities completely different from 
those reported by Sofue & Yoshida (1993). Hence, our de- 
tection is more in agreement with the former authors, given 
their level of sensitivity. 



3 DISCUSSION 

3.1 Characteristics of the CO emission 

The spectra displayed in Fig. ^| show that the 12 CO (1-0) and 
12 CO(2-l) emission is detected in the heliocentric velocity 
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Table 1. Line parameters for the complex: (v) is the heliocen- 
tric centroid velocity and o"fwhm gives the measured FWHM of 
velocity profile. The CO flux, Sco, an t the CO luminosity, Iqo 
(expressed in units of T^ , the antenna temperature corrected for 
atmospheric absorption and rear sidelobes) are given for both 
lines, rja is the the aperture efficiency. 

12 CO (v) o-fwhm -fco Sco Va 

(km s" 1 ) (K km s" 1 ) (Jy beam -1 
km s -1 ) 



(1-0) 
(2-1) 



-156 
-154 



-35 
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0.36±0.05 
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Figure 4. Extinction curve in the dust complex. Each symbol 
corresponds to an average over 854 pixels, with an extinction 
measured in BVRI and for which A/ > 0.04. They show the 
extinction law estimated for different geometries of the dust. The 
lines correspond to the models of Fitzpatrick (1999). 



3.2 The extinction curve for D395A/393/394 

The extinction curve is determined by comparing the ex- 
tinction in B, V , R and J on a pixel basis. Figure^ shows 
that the Av/E(B — V) ratio Ry- is closer to 2.0 than to 
the canonical value of 3.1. This cannot be accounted by the 
geometry of the dust distribution (see Walterbos & Kenni- 
cutt (1988) for a comprehensive review). If dust and stars 
are uniformly mixed along the line of sight, the true ex- 
tinction would then be increased by a factor of 2, but the 
extinction curve would still present the same Rv. If the 
dust is located in a thin layer, which is rather plausible 
if associated with the molecular and ionised gases, the ob- 
served extinction, as defined in Sect, 
to — 2.5 log 10 (:r + (1 — x) x exp(— r)), where 



then corresponds 
t is the real 

optical depth at a given wavelength and x the fraction of 
light in front of the dust. As observed extinctions as large as 
0.49 mag are measured in B, the fraction of foreground light 
x has to be smaller than 0.6. This fraction is 0.3 and 0.02 
if the dust patch lies in a plane inclined by 45° and 77° re- 
spectively, when using Galactic bulge and disc models from 
Binney, Gerhard & Spergel (1997). Similar deviations from 
the average Rv have already been reported in a few lines of 
sight across M31's disc (Massey et al. 1995, and references 
therein): like what is observed in the Galaxy (e.g. Szomoru 
& Guhathakurta 1999), local fluctuations are also present in 
M31. 

For x = 0, we estimate the extinction Ay=0.06 mag 
beam -1 when averaged over the 115 GHz beam, correspond- 
ing to E(B — V)=0.0285 mag beam -1 . Given the arguments 
discussed above, the uncertainties affecting these values are 
negligible for the following. 



range Vkelio = —130, —190 km s . The measured centroid 
velocities are in good agreement with the velocities of the 
ionised gas published by Boulesteix et al. (1987). This sug- 
gests that the molecular clouds lie in the same plane as the 
ionised gas; the latter is thought to be inclined by 45° in 
the region of interest (see the discussion by Ciardullo et al. 
1988). At the position observed, the ionised gas presents a 
velocity gradient of ~30 km s -1 over the 2l"-wide telescope 
beam. This is consistent with the FWHM of the CO (1-0) 
line profile (31 km s -1 ) measured for the molecular com- 
plex. This suggests that we are observing not one single, 
but several clouds distributed across the telescope beam. 

The comparison of the intensities expressed as the 
12 CO(2-l) to (1-0) line ratio R is not straightforward as the 
beams are different and the geometry of the emitting regions 
is a priori unknown. If the emitting region were uniform and 
much larger than the beams, the line intensity ratio would 
simply be R=T^(2-l)/T^(l-0)=0.6. We adopt the more 
realistic assumption suggested above that the dust follows 
the gas distribution. Integrating the extinction distribution 
(Fig. ^ over the telescope beam and taking into account the 
difference in telescope beam shapes and efficiencies between 
115 and 230 GHz (Greve, Kramer & Wild 1998), we arrive 
at a ratio R=0.5. This value is typical of those observed in 
M31's molecular arms (Neininger et al. 1998) but is larger 
than those observed by Loinard et al. (1996). 



3.3 The gas-dust connection 

Although it is very dubious to interpret the observed CO 
luminosity in terms of mass of molecular hydrogen, the 
Ico /E(B — V)=9.8 K km s -1 mag -1 ratio measured here 
is remarkably similar to the value observed for the molec- 
ular clouds in our Galaxy. Hence, the most direct method 
to get the order of magnitude for the mass of the gas in 
this complex is to assume a standard dust-to-gas mass ra- 
tio, N H /E(B - V) = 5.8 10 21 atoms cm -2 mag -1 (Bohlin, 
Savage & Drake, 1978). Integrating directly over the distri- 
bution of extinction, we derive a mass of 1.5 10 4 M (10 4 M 
in the main solid angle of the beam at 115 GHz). This mass 
could be overestimated if the metallicity is larger than 1 in 
the central region of M31. Continuing, based on the optical 
extinction map, we find that the size of this complex is about 
60 pc; this is typical for a giant molecular complex. However, 
its average volume density (~ 1 cm -3 ) is extremely low: two 
orders of magnitude lower than the typical densities in gi- 
ant molecular clouds. For this reason, this complex must be 
highly clumped as supported by the detection of these CO 
lines. 

If these clumps are gravitationally bound together, we 
estimate a 
virial mass Mvir/M© ~ 200r/(pc)((JFWHM/(km/s)) 2 , i.e. 
Mvir ~ 1.6 10 7 Mq. The difference of three orders of mag- 
nitude compared to the above result based on the optical 
extinction is so large that either these clumps are not bound 
together or the mass is not dominated by the gaseous com- 
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ponent. For the very same reason, this complex does not 
follow the velocity line-width versus diameter relationship 
found by Solomon et al. (1987) for Galactic molecular clouds 
near virial equilibrium, but presents similar characteristics 
to the clouds studied by Oka et al. (1998) in the Galactic 
centre. 



4 CONCLUSIONS 

Emission of the (1-0) and (2-1) lines of the 12 CO molecule 
has been detected near the centre of M 31 at 15a and 10<r 
respectively. The CO centroid velocity of —155 km s~ sug- 
gests that this molecular gas could be located in the ionised 
gas disc detected in these regions. The observed velocity 
dispersion of the molecular gas in the beam is compatible 
with the velocity gradient of the ionised gas. Together with 
the patchy appearance of the extinction map, this indicates 
the presence of several clouds distributed over the beam. 
The lco/E(B — V) ratio is remarkably similar to the val- 
ues observed in molecular clouds in the Galaxy and M31's 
disc. Based on a standard gas-to-dust ratio, the mass of the 
molecular complex is of order 10 4 Mq. In order to reach 
densities compatible with the excitation of the CO lines, 
the complex must be highly clumped. As previous studies 
with poorer resolution have shown, the mass inferred is small 
compared to the mass content of the Galactic centre (Mor- 
ris & Serabyn 1996). With such a mass the line- width of 
the CO line is by far too broad for a gravitationally bound 
complex, which explains previous inconsistencies based on 
this assumption. 

Accordingly, we modelled the extinction with a random 
distribution of small spherical clumps, all identical with a 
size r and density p. They all lie in a sphere of radius 13", 
tracing a total mass of 10 4 M Q . Models with ~250 clumps 
reproduce the main features of the histogram of measured 
extinctions. With the previous assumption of a common lo- 
cation for the molecular, ionised gas and the dust, the con- 
figuration (r,p) which reproduces best the measured extinc- 
tion correspond to clumps with r~2.6" and p ~400 H cm -3 . 
Following the LVG homogeneous cloud models of Garcfa- 
Burillo, Guelin & Cernicharo (1993), the measured line ratio 
corresponds to CO (1-0) and (2-1) line excitation tempera- 
tures and to a kinetic temperature all >10 K. This modelling 
will be further discussed in a subsequent paper. 

The "standard" CO-luminosity to visual extinction ra- 
tio and the relatively high CO (1-0) excitation temperature 
seem to rule out the presence of large amounts of hidden CO 
and H2 in the inner bulge of M 31, in apparent disagreement 
with the result of Loinard et al. (1996) in the inner disc. 
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